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‘What happened to the Dar sping? 
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Itis well known that you can measure the 
Synchronous Ampliication Factor of a rotal- 
lng machine at some reasonable value, like 4 
to, and then have the machine get into a bad 
forward circular Instability whip at operating 
speed and conditions. Por such an instability 
to occur, the Operating Amplification Factor 
must have gone to infinity. 

The obvious question ts: “What happened 
to the Damping?” k will be shown here that 
the damping 1s in good shape (more so than 
you thought) and that it was the Swirling 
partion of the Quadrature Stiffness that ate 
your lunch, 

Consider a rotating: machine that operates 
above at least one of its scll-balance reso- 
nance speeds (“criticals’). In order to evalu- 
ate its net effective Spring, Damping, Mass, 
and its Synchronous Amplification Factor, 
the usual process 1s to first get the rotor well 
balanced. Then add an imbalance of a partic- 
ular weight, radius, and angle at a lateral 
position that will excite the balance reso- 

Next, run the machine from zero speed to 
a speed well above the balance resonance 
speed. Document the rpm, amplitude, and 
phase of the rotative (synchronous) speed on 
start-up or coastdown across the specd 
range. Subtract the Slow Roll bow, and plot 
the data in Polar and Bode plots. 

For this test, the machine should be as 
nearly as possible in its operating configura- 
tion, load, and environment. This applies to 
its foundation and alignment as well as to the 
characteristics of its energy source (in the 
case ofa driver element like a gas turbine), or 
the characteristics of its load fin the case of a 
pump, compressor, and generator). 

Now, let's work an actual problem for a 
typical rotor system, Here are the pertinent 
characteristics of a typical machine (except 
one vital parameter, which wall be shown 
later). The net spring coefficient is 1 million 
lbs./in. The weight of the rotor is 3300 lbs., 
and its effective weight (fit were all concen- 
trated in one spot at midspan) is 2316 Ibs. 
Dividing by gravitational acceleration of 386 


in.sec,* , the rotor effective mass Mupp 18 6 
lbs. sec. +n. 


The operating speed is 7640 rpm (which, 
to make calculations easier, comes out to #00 
rad sec.j. From this information, the first 
balance resonance speed can be found as 
follows: 


Opes = VK = W 10% = 408 radusec. 
Mier 6 


or 


RPMpes = ( Oy Wpps = 3900 rpm 
fay 


Figure | shows the rotor configuration, its 
lirst mode bow shape, and the location of the 
shalt observing probe. 

By looking at the Polar plot or the Slow 
Roll Compensated Bodé plot (Figure 2), you 
can work out the Synchronous Amplification 
Factor, and Observed Damping 


she pole pee —_ wpro X Mew 
QeyNcH 


If you like to use synchronous damping 
coefficient zeta, 1 1s: 


Eu YNCH — l 
Wyne [ 


and, for logarithmic decrement buffs, it is 
the 


Synchronous Log Dee = 
QeyNCH 


The reason for the subscript adjectives, 
according to my colleague Joe, is that a 
rivaling machine can do more tricks than a 
monkey on a 100 foot (30.48 meter) rope, 
Therefore, don’t throw away the subscripts. 

From Figure 2, you can figure out the 
Synchronous Amplification Factor by the 
half power method, the rato method, or the 
phase lag rate method. All these methods 
work okay, but strictly speaking, you should 
divide amplitude by speed before doing the 
half power method. 


AS you can see, the Qesyrgepy 15 about 7.0. 
From this you can get 
Synchronous Damping Ratio 
e=_! = 004 
aQ 


E 


Synchronous Log Dec = 0.449, and 
Observed Damping 


— (408016) — 450 Ih coe 
Dons = s = 350 Ib secän 


All finished? Not by a long shot. You next 
need a new graph. Nobody is using this 
eraph yel—cacepl my colleagucs— Joe and 
Agnes — and me. We think it is extremely 
Unportant because it gives the engineer a 
clear picture of rotating machinery behavior. 
Come to think of it, the rotor stability 
researchers al Cal Tech and Tokyo University 
also use this seraph. 

Anyway, this is pretty simple mechanical 
engineering, To make this graph, you simply 
write out the total of all the dynamic stiffness 
terms that you have so lar: 


Observed Dynamic Stitiness, 


Koss =K- 


w? Mere + jwDons 

This is a vector equation; (1) because | pul 
a bar on top of Kags, and (2) because, 
whereas, K and —ow@M take off at 0° and 
180°, respectively, the + jwD,,),. 15 at righi 
angles (which are much better than lett 
angles, or wrong angles, according to Joc. | 
am lett handed, so 1 sorl of resent that. 

Instead of fooling with this vector, break it 
up into two parts, the Direct happenings at 9° 
and 180°, and the Quadrature happenings al 
+90". Thus, there are two graphs, Kray eer 
and Kouanraruer Loe has a bg problem 
with the word quadrature: he thinks it has 
something to do with the number 4. He's tn 
good shape, however, compared to some 
people who call any term with a W | 
imaginary!) 


Thus, Koes 


CO ADRATURE: 


Kons DIRECT + JEons 


DE -i = i aes 2 
and Koss precr = K 7 w*Mgrr 
and kons OUADRATURE ` WD ops 


The values are: 
K = 10 Ib.fin. 
Dong = 330 Ib. secin, 
Mefr = 


wis Oto 1000 rad. sec. (0 to 9550 rpm) 


6 lb. sec. “fin. 
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The graphs of the Direct and Quadrature 
Dynamic Stiffness are shown in Figure 3, 

There are more terms in this equation for a 
typical rotor system, some of them linear, 
some nonlinear, some symmetric, some 
nansymimetric, some having to do with 
mechanical things, and some to do with 
aerodynamics or hydrodynamics. 


The list of these is impressive, including: 
UPPOO ASY TImETY 

* Kotor asyminetry 

* Fluidic inertia 

e CYTOSCOpLcs 

© Swirling rat of liquids and gasses 

© Pull and partial rubs 

‘Coupling lockup 

* Lateral effects by torsional action 

” Shifts of bearing loading due to alignment 
* Shilts of bearing loading due to load. 


However, even though these and other 
terms might be influencing the Dynamic 
Siffness of a rotor system, the swirling ratio 
af the liquids and gasses is clearly the big one 
to watch. 

The swirling ratio, A (lambda), occurs on 
nearly all machines. [tis the net ratio of the 
rotation rale of steam, liquid, andor gasses 
inthe machine around any and all “cylinders 
inside cylinders.” Therefore, bearings 
(except tilting pad and rolling element bear- 
ings), virtually all seals. steam shrouds, all 
parts of pumps, and onl retainer rings contrib- 
ute to swirling. 


Damping Stiffness is a Dynamic Quad- 
rature term because its term 1s viscous damp- 
ing times velocity, and velocity is 90° ahead 
of displacement motion. The Swirling 
Dynamic Stiffness Is also a Quadrature term, 
but for an entirely different reason. The 
Swirling Dynamic Suliness is the product of 
three terms: (1) the net swirling ratio, (2) 
rotative speed, and (3) viscous damping. 
Swirling liquids and gasses generate an 
active force on the shalt when they are 
trapped by a wedge. The wedge is sel up 
behind the constriction formed by the side- 
wise shaft motion, Because of this, the 
Swirling Dynamic Stiffness has the form 


Kwp =L E (ANG perp?) * iD) lb.sin. 


This term is generally referred to as a 
“cross spring, or" KY. or Kyy oraera- 
dynamic cross coupling. A better term, 
when dealing with sleeve bearings, is “oll 
wedge support stiffness.” 

The dynamic stiffness of viscous damping 
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is (+ [Ors epee? (Dh which is the shaft pas- 
sively pushing on the fluids or gasses. The 
dynamic stiffness of the swirling 1s 
(—[M@pory! * (2), which is the fluids or 
gasses actively pushing on the shaft. 

These two terms are the principal contrib- 
utors to the total Quadrature Dynamic Stiff- 
ness, so that 


Kouan = {— Mape YEN + Wp pre D) lb.ín, 


Therefore, one more parameter is needed 
to complete the description of this rotor sys- 
tem, the average swirling ratio, A. A very 
gowl assumption for this machine, and for 
that matter, any turbine, compressor, pump, 
or generator is about 0.40 to 0.45. (Joe 
thinks that this is related to regular old oil 
whirl, and he is cuht!), Occasionally, rotors 
whirl al that frequency. 

With this new pece of information, Fig- 
ure 4 shows the actual behavior of this 
machine, assuming that A = 0.41. Note that 
the horizontal axis is labeled “Precession 
Speed Wpppc- [he Direct Dynamic Stiff- 
ness is unchanged because esyroscopics and 
fluidic inertia, which contribute to the Direct 
Stiffness, are not immediately important, but 
there is a tolally new Quadrature Dynamic 
Stillness term. 

The old rule tor figuring damping, 
(png Mpg) 

synca? 
does net give the correct result, unless the 
swirling is zero! 


D= 


Here is what works with excellent 
accuracy on rotating machinery: 
— (pps Mere) 
(QeywcH la) 


_ _ (4086) _ 393 Ib. sesin. for 
(Fa 1-0.41) this machine 
example. 


The viscous damping is absolutely NOT 
the 350 Ib. sec.‘in, previously calculated! 
The good news ts that the damping 1s much 
higher than that measured by the data from 
start-up or rundown. 

The bad news is that when the rotor is at 
7640 rpm (800 rad./sec.), the Swirling Stiff- 
TESS Is: 
Kwp = (O41 80593) = -195,000 
Ib. fin, 

As a result, the Quadrature Dynamic Stiff- 
ness at balance resonance (first critical) 
when this machine 1s at tts operating speed 
of 7640 rpm (800 rad. sec.) ts: 
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Kew Œ 3900 rpm 
-195 OOK) + (S95 (408) 
-195 (KK) + 342 000 
+ 47 000 Ibin. 
instead of the 143,000 |b./in, 
observed on start-up, 

The Operating Amplification Factor is: 
(w Res M ggr) 


that was 


Qop = 
Kouan 
_ (408 2x (6 _ 9 
A ON) 
When the machine 1s at its operating speed 
above resonance, this is the vital data neces- 


sary, nol the Qeyeageqy OFT that appeared on 
the start-up data. 

This behavior is very typical of all rotating 
machinery. Figure $ shows the amplitude 
and phase response as the machine sees its 
dynamic impedance when it is al operating 
speed. Compare this to the amplitude and 
phase observed synchronously in Figure 2 


(dashed line in Figure $5) and the necessity of 


the Direct and Quadrature Dynamic Stiff- 
ness plots becomes obvious. 

The necessity of specifying the Operating 
Amplification Factor is also obvious. The 
example used here has a very slim marein of 
stability. (When the Direct and Quadrature 
Dynamic stiffness terms are both at 0 at the 
same precession speed, the machine will be 
whipping at that speed.) With these graphs, it 
is simple to calculate how much loss of 
spring stiffness K will introduce instability 
(Joe says that if this machine goes from | x 
10° Ibvin. to K = (O.404¢8000-(6) = 
645,000 [bs in., the machine will be in 
whip. } 

Önce this knowledge and method comes 
Into general use, the next step will be to take 


up the higher engineering level problem of 


controlling the machine parameters to pre- 
cictably chimmate instability, Some immedi- 
ate clear cures of instability include: (1) 
design out the swirling, and (2) if you can't 
do that, hold the roter at high eccentricity 
ratio, like O.7 or so, in seals and cylindrical 
bearings with steady-state loading. 

For copies of stability papers Dr. Agnes 
Muszynska and [ have published, please 
contact our Customer Information Center, 
P.O. Box 157, Minden, NY 89423, (702) 
782-3611. E 
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Figure 2 


_ CES ERED MECT STP RRS 
| 
| 
l 

g ] 
i 
$ F 
"| 
HE) za] | 
l 
HGG i 





| 
eagoa000 ee 
| 
ri Hi fae i 
Ci) peng radan 
SAFT ROTATIVE SPEED 
=i, 


Figure 3 
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Figure 4 
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